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I. INTRODUCTION 



(Dated: February 4, 2008) 

Quasinormal modes for scalar field perturbations of a Schwarzschild-de Sitter (SdS) black hole 
are investigated. An analytical approximation is proposed for the problem. The quasinormal modes 
are evaluated for this approximate model in the limit when black hole mass is much smaller than 
the radius of curvature of the spacetime. The model mirrors some striking features observed in 
numerical studies of time behaviour of scalar perturbations of the SdS black hole. In particular, 
it shows the presence of two sets of modes relevant at two different time scales, proportional to 
the surface gravities of the black hole and cosmological horizons respectively. These quasinormal 
modes are not complete - another feature observed in the numerical studies. Refinements of this 
model to yield more accurate quantitative agreement with numerical results are discussed. Further 
investigations of this model are outlined, which would provide a valuable insight into time behaviour 
^"v of perturbations in the SdS spacetime. 
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A characteristic feature of the response of a black hole to external perturbations is the appearance of quasinormal 
modes (QNMs). The presence of these modes was first noted in a study of perturbations of the Schwarzschild black 
hole, by Vishveshvara p|. Since then, QNMs for asymptotically flat black holes have been computed by a variety of 
(f-) ' numerical and analytical approximation methods. A detailed review of this work can be found in 

o: 

""q I QNMs were first found in a stability analysis of black holes. However, the observation that these damped oscillations 
Q" 1 are intrinsic characteristics of the black hole exterior geometry and depend only on the black hole parameters implies 
that they are the imprint of a black hole in its response to perturbations. Further, QNMs are even seen at intermedi- 
(3Tj, ate/late times in fully non-linear situations like systems undergoing gravitational collapse. Thus they are expected to 
play a significant role in the search for gravitational waves and black holes. Recent evidence for a non-zero positive 
cosmological constant points to the importance of studying black holes in such a background. The simplest black hole 
in this class is a Schwarzschild black hole in de Sitter space (SdS). The QNM spectra for gravitational perturbations 
of this black hole have been investigated in [J, |5j using numerical and analytical approximation techniques. QNM 
spectra for the SdS black hole in the near-extremal case when the two horizons are nearly coincident (i.e, the black 
hole mass and radius of curvature of the spacetime of the same order) are derived in |(| . An interesting numerical 
study of evolution of scalar fields in the SdS spacetime has been performed in 0, @ • This numerical study reveals 
the differences between the response of the Schwarzschild black hole and the SdS black hole to perturbations in the 
physically interesting regime where the black hole mass is much smaller than the radius of curvature of the spacetime 
(which is a relevant regime because the cosmological constant is small). The SdS black hole shows QNM behaviour 
similar to the Schwarzschild black hole at intermediate times, and soon changes to a power law decay. At late times 
it again shows QNM behaviour (exponential decay) where the QNMs are now proportional to the surface gravity of 
the cosmological horizon. For the Schwarzschild black hole, analytic approximations have been used to give a general 
(approximate) formula for the QNMs as a function of the black hole parameters and the angular mode 0- It is of 
interest to compute, in the same spirit, the QNMs for the SdS black hole in an analytical approximation scheme in 
the physically relevant regime. The advantage is that there exist known numerical studies which reveal characteristic 
features in the time decay of a field in the SdS black hole. The analytical approximation could be refined to reproduce 
these features, and at the same time would give a general expression for the QNMs as a function of the spacetime 
parameters and the angular mode. The approximate model, being simpler to study may also be possible to analyse 
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completely. It would thus give insight into the qualitative features of the complete time behaviour of the field in SdS 
spacetimc. Further, the near-extremal limit and its effect on the QNMs and time behaviour of fields could be explored 
in this model. 



In the recent past, the QNMs of Anti-de Sitter (AdS) black holes have been studied extensively. A detailed numerical 
study of QNM decay of scalar fields in AdS black hole backgrounds in various dimensions was performed in 01 • Exact 
computations of QNMs for the BTZ black hole in (2 + 1) dimensions were demonstrated first in ^l| and subsequently 
higher order modes (as also numerical results for gravitational perturbations of SAdS black holes) are shown in (12j . 
More general numerical computations of QNMs for higher dimensional AdS black holes are found in [13. Decay of 
scalar fields coupled to curvature in topological AdS black hole backgrounds is studied numerically in |l4. flU and 
analytically in (lfi|. An interpretation of the QNMs for an AdS black hole in the dual CFT using the AdS/CFT 
correspondence is provided in [13 (see also [18|]). An interesting question is whether the QNMs of the SdS black hole 
have an interpretation in a dual CFT using the recently proposed dS/CFT correspondence 0] (for a study of pure 
de Sitter QNMs in this context, see [2(| ). The study of the time behaviour of perturbations in an analytical model 
for the SdS spacetime could be used to study such a possible interpretation. 



In this paper, we propose an analytical model to study the time behaviour of perturbations in the SdS black hole. 
We compute the QNMs of this model in the regime when the black hole mass is small compared to the radius of 
curvature of the spacetime. The SdS spacetime is bounded by the black hole and cosmological horizon. Our analytical 
approximation consists of approximating the SdS potential by a potential that reproduces its asymptotic behaviour 
at both boundaries. In section 2, we enumerate the QNM boundary conditions for the SdS black hole, and describe 
our approximation potential along with the motivations for our choice. The chosen potential has a discontinuity in 
its second derivative at the maximum. In section 3, we describe the exact solutions to this potential and list the 
correct matching coefficients at the maximum of the potential. The actual computation of the QNMs is done in the 
regime of interest in the next section. Two sets of QNMs, one proportional to the black hole horizon surface gravity 
and the other proportional to the surface gravity of the cosmological horizon are found - the first set describing the 
time behaviour of the field at certain intermediate times, and the latter describing late-time behaviour. A qualitative 
picture for the time scales where each of these two sets of QNMs comes into play is given - and it is discussed why 
this behaviour is expected only in the regime where the black hole mass is small compared to radius of curvature. 
The real part of the QNMs depends on the maximum of the potential. Our approximation consists of identifying 
the value of the maximum of the model potential with that of the SdS potential. We compute the maximum value 
of the SdS potential in section 5. This unfortunately is difficult to do for any general angular mode. We therefore 
consider the small angular mode and large angular mode cases separately while computing the SdS maximum. In 
our last section, we present a detailed comparison of the results of this model with numerical studies. There is good 
qualitative agreement in the most striking feature of the numerical studies, the presence of the two sets of QNMs 
relevant at two different time scales, and which are proportional to the black hole and cosmological horizon surface 
gravities respectively. We also show in this section, non-completeness of QNMs in our model - that is also expected 
from the observation of power law behaviour at intermediate times for the SdS black hole. We discuss how the model 
could be refined or modified to reproduce the more quantitative features of the numerical results like the dependence 
of the QNMs on the angular mode. We also discuss the exciting questions which can be addressed for this model, and 
would shed light on our understanding of perturbations of the SdS spacetime. 



II. ANALYTICAL APPROXIMATION OF THE SDS POTENTIAL 



The metric for the Schwarzschild-de Sitter spacetime in four dimensions for cosmological constant A = 3 /I 2 is given 

by 

ds 2 = -{Nf dt 2 + (N)- 2 dr 2 + r 2 (dfl) 2 , (1) 

with 



2M r 2 N 

iV = V (1 - — -F J - (2) 

There are now two horizons given by the zeroes of N 2 . These are the two real positive roots of iV 2 , 

2 1 

?'b = —= sin[- sin _1 (3AfV / A)] and 
VA 3 
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r c = ^sin^shr^MVA) + y] . (3) 

The third root of N 2 is negative, and given by ro = — (r& + r c ). rb is the black hole horizon and r c the cosmological 
horizon. As in the asymptotically flat case, the spacetime has a curvature singularity at r = 0. 

We consider a massless scalar field $ in this background. Then, the Klein-Gordon equation for the field is 

d^V^gg^h)* = o . (4) 

Using an ansatz for the field 

$ = -x(r)e i " t Y LM (9,<t>), (5) 
r 

and by going to the tortoise coordinate x where dx = dr/N 2 , the Klein-Gordon equation is 



pt + V L (r) X =u J 2 X , (6) 



dx 2 
where 

When r = r/,, a; = — oo. When r = r c , x = oo. We are interested in studying scalar field perturbations in the region 
bounded by these two horizons. The potential Vl(t) goes exponentially to zero as a function of the x coordinate as 
one approaches the two horizons. However, the rate of approach is not the same - it depends on the surface gravities 
associated with the horizons. The surface gravity at a horizon is defined as a/, = | \df/dr\ r . Thus the surface 
gravities associated with the black hole and cosmological horizons, ctb and a c respectively are 



(r c - r b )(r b - r ) 

2l 2 r b 
(r c - n)(r c - To) 



a b = and (8) 



2l 2 r r 



(9) 



The potential Vt(r) goes to zero as e 2abX as x — > — oo, i.e as the black hole horizon is approached. When x — * oo, 
i.e as the cosmological horizon is approached, Vl(v) goes to zero as e~ 2oicX . In studying the time evolution of the 
perturbation, it is important to first compute the QNMs for the problem. QNMs are given in this case by solutions to 
© with specific boundary conditions - ingoing at the black hole horizon, and outgoing at the cosmological horizon. 
These are the physically motivated boundary conditions for the problem as a timelike observer in this spacetime can 
send messages, but not receive them through the (future) cosmological horizon. The boundary condition at the black 
hole horizon is the usual one for computing QNMs of black holes and represents the classical absorption by the black 
hole. The QNM boundary conditions for the field are thus: 

X oc exp(zo;x); x — * — oo and 

X oc exp(— iux); x — > +oo . (10) 
lo is now complex and lj — loq + ir, where T > 0. 

Thus, the QNMs can be obtained by solving JfJJl with the above boundary conditions. However, exact solutions 
to this equation are not known. In fact, such exact solutions are not available even for the Schwarzschild black 
hole in an asymptotically flat spacetime. Instead, what is done in the asymptotically flat case is to either solve 
the equation numerically and compute the QNMs, or to analytically approximate the Schwarzschild potential by an 
exactly solvable potential. This potential is described by parameters which can be used to fit the potential to the 
Schwarzschild potential . We follow the second approach of an analytic approximation in computing the QNMs of 
the SdS black hole. 
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In our search for a potential that would be a good analytic approximation to the SdS potential, we set out to 
reproduce the behaviour of the SdS potential as x — > ±00. The approximate potential V that we have chosen has the 
following features : 

1. Its second derivative is discontinuous at r = r m , where r m is the value of r for which the actual SdS potential 
has a maximum, i.e a solution of cLVl / dr = 0. 



2. In the region r b < r < r m , the potential is 



[cosha b (x - x m )] 2 

x m is the value of the x coordinate when r = r m . is given by ©■ V\ is referred to as the Poschl- Teller potential 
in literature, and was first used in [t| to approximate the Schwarzschild potential. V m is taken to be the maximum of 
the SdS potential. 



3. In the region r m < r < r c , the potential is taken to be again the Poschl- Teller potential, but with a different 
choice of parameter. It is 



V 2 = . (12) 

[cosh a c (x - Xm)] 2 



a c is given by ©. 



As can be checked, this potential reproduces the asymptotic form of the SdS potential at the two boundaries. 
Furthermore, this potential is exactly solvable. Therefore what remains is a computation of V m , the value of the SdS 
potential at its maximum. Computing the maximum of the SdS potential is difficult as the potential is complicated. 
We are forced to work in the regime M < 1. The Poschl- Teller has been used to approximate the SdS potential 
recently [f|. However, in this case, one Poschl- Teller potential is used throughout - and this is a good approximation 
only for the nearly extreme SdS black hole, where the two horizons are very close (i.e, M and I are of the same order) 
0. In the other regime M -C I, one Poschl- Teller potential provides a poor approximation to the SdS potential as the 
asymptotic forms of the SdS potential at both the boundaries cannot be reproduced by it. This problem is solved by 
our choice for the potential, V. It would be exciting if the near-extremal limit in our model led to the QNMs found 
in HJ. We must now check what features, qualitative or quantitative, we reproduce of numerical studies of dynamics 
of a scalar field in the SdS spacetime done in 0,0 by a study of this potential. 

We now proceed to compute the QNMs of the approximate potential V which is given by V\ for r < r m and V2 for 
r > r m . We are interested in the solutions to the equation : 

" S + VX = ^ X ■ (13) 

For both V\ and V% , exact solutions can be found. We first examine the exact solutions to I1K.'>I> with both V = Vi and 
V = V2. Then, the solutions have to be matched at r = r m . Finally, we must pick the solutions that obey the QNM 
boundary conditions at both the boundaries. 



III. EXACT SOLUTIONS FOR THE POTENTIAL 



Let us consider the solutions to with V = V\. The solutions to this (Poschl- Teller) potential are known. They 
can be written in terms of the hypergeometric functions after a change of variables to £ where £/(l — £) = e 2a b(x-x m ) _ 
Then, two linearly independent solutions are : 

Xi = [£(l-£)]^F(a,6,c;£) and (14) 

X2 = [tfi-O^CO^fo-c+M-c+M-cO. ( 15 ) 
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We have 




1 1 V m iu vjj 

o~ 1/7 2+ — ■■ c= H 

2 y 4 at a b a b 



2 • ■ - 1 + — ■ ( 16 ) 



The black hole horizon is at £ = 0. At the maximum of the potential, when x = x m , £ = 1/2. 
Solutions to (|13|) with V = V2 are similar to the previous case, but with a b replaced by a c . 

Define fi where £i/l - £1 = e 2a ^ x ^ . 

Xi = Ki(l-6)]^^(«i,fei,ci;ei) and (17) 
X2 = [6(1 -6)]^(a) 1 " Cl ^(ai-ci + -a + 1,2 - Cl ;6) ■ (18) 

Here 

1 . /l V m iw 1 1 V m iuj iuj 

ai = - + W- 5- H , 61 = - - W - 5- H , ci = H . (19) 

2 y 4 a 2 a c 2 y 4 a 2 a c a c 

The next step is to match an arbitrary superposition of the two solutions l|15|) and l|15fl with an arbitrary superposition 
of (|18H and i|18|) - at x = x m . This involves matching both the wave functions and their first derivatives across 
x = x m . Let us use the following notation: Xi( x m) = A, Xi(x m ) = Y. Also, ^-(x m ) = A', ^r(x m ) = Y'. 
Xi(x m )=X, x 2 (x m )=Y : and ^(x m ) = X', *&(x m )=Y'. 



Then matching the wave functions and their derivatives across x — x m yields 

AX + BY = AX + BY ; AX' + BY' = AX' + BY' . (20) 

Here, A, B, A and B are arbitrary coefficients of superposition of the two linearly independent solutions, which are 
then x — A\i + B\2 for r < r m and \ — Axi + £>X2 for r > r m . 

A and B are related to A and B as : 



r , (XY'-X'Y) (YY'-Y'Y) , . 

A = A^= ^-J- + B -W^ , (21) 

(XY'-X'Y) (XY'-X'Y) 

, (xx' - x'x) „ (rl' - r'l) 

(YX'-Y'X) (YX'-Y'X) 

We recognise that at x = x m , the value of £ that appears in (|15|l and (|15|l is 1/2. Using this and properties of the 
hypergeometric functions, we can write : 

A = V^Fexp ) — 2 , 2 , , 23) 

V( a b ^r(i + f)r(i + |) 

X' = A^exp(-^)^ F(2 + ^ . , (24) 

2 PV « fc ; c r(i + f)r(i + |) v ' 

_ iuln2 r (l-^) / ok n 
F = V7rexp( ) 7— , (25) 

v PV a b ; r(i-f)r(i-|) ' v ' 

v i ■ / — ,i^ln2 r (l^^) 

+ f (a ~ C+ 9 1)(b " C + 1) F(a-c + 2,6-c + 2,3-c;l) . (26) 
2 2 — c 2 
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X, Y, X' and Y' are given by the same expressions as above, except that we replace aj, with a c . Consequently a, b 
and c are also replaced by their counterparts a%, b\ and c\. We therefore do not display those expressions here. 



IV. COMPUTATION OF QUASINORMAL MODES 



We search for solutions to the problem satisfying quasinormal mode boundary conditions at the black hole and 
cosmological horizons given by (|10[) . The solution obeying the QNM boundary conditions at the black hole horizon 
Xqnm for x > x m is given by A\i + B\2- A an d B for Xqnm are given by i|21|) and 1|22|) with B = 0. Near the 
cosmological horizon, i.e as x —* oo, this solution is 

Xqnm ~ AP\ exp(— ibjx) + AQ\ exp(iu>x) + BP2 exp(—iuix) + BQ2 exp(iu>x) , (27) 

where 

p _ r(ci)r(ci - oi - 61) 

I (ci - ai)l (ci - 61) 
= r ' C -' r '-. r + ^.- C - ) , (29, 

r(ai)r(oi) 
_ r(2- Cl )r(ci-ai-bi) 

2 — — p?i vp?i — ' 

r(i - oi)r(i - hi) 
n _ r(2- Cl )r(a 1 + b 1 -c 1 ) 

^ 2 r(ai - d + i)r(6i - Cl + 1) ' 1 j 

For l|27|) to obey QNM boundary conditions at the cosmological horizon given by (|10f> . we must have 

igi + BQ 2 = 0, (32) 



Substituting for A and -B from l|21ll and (|22|l with _B = 0, the above equation is 

(XY-X'Y) T ^ (XX'-X'X)- ^^-^ —— = . (33) 

r(oi)r(6i) r(ai - ci + i)r(6i - ci + 1) 

The special frequencies cj that solve the above equation are the quasinormal mode frequencies for this problem. 
However, as can be seen from the expressions for X and X' from lIl'MI) and l|24l) - the zeroes of X and X' do not occur 
for the same set of frequencies. The zeroes of the above equation cannot therefore be computed simply. However, here 
we reflect on what we must expect. From the numerical work in 0, Q , we note the presence of two different patterns 
of exponential decay of the field with time. For intermediate times, the field decays exponentially with the exponent 
proportional to the black hole surface gravity. This is a typical QNM time behaviour where the QNM frequency is 
proportional to the black hole surface gravity. This temporal behaviour is similar to that of a field in a Schwarzschild 
black hole background. This exponential decay is followed by a power law decay of the field - but for late times, there 
is a switch back to exponential decay. Now the exponent is proportional to the surface gravity of the cosmological 
horizon. This exponential decay is again a QNM time behaviour, but the QNMs are now proportional to the surface 
gravity of the cosmological horizon. This numerical result was seen for a choice of well separated length scales M and 



We expect a similar behaviour (i.e two sets of QNMs) for our approximate potential V . At very early times, there 
is a direct propagation of the perturbation without scattering by the potential. At later times, the perturbation gets 
multiply scattered by the potential. It therefore starts reflecting the properties of the potential. The time behaviour 
of the field here is dominated by the QNMs. In our case, the second derivative of the potential has a discontinuity 
at r = r m . The QNMs which we have obtained (and which are supposed to describe the time evolution after the 
very early times) also reflect this. When the perturbation gets scattered for rj, < r < r m by the potential (which is 
now Vi), its time evolution will be governed by QNMs related to V% which should be proportional to ctb (as is typical 
for QNMs of a Poschl Teller potential). A rough estimate for the times when this occurs is the light crossing time 
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to traverse a distance of the order of r m , i.e r m /c. When the perturbation travels a distance much greater than r m , 
it starts getting scattered by the potential V2, and its time evolution is governed by QNMs related to V2, which we 
expect to be proportional to a c . This roughly occurs at times of the order of the light crossing time to traverse a 
distance I, i.e l/c. At such times and later, the decay of the perturbation is governed solely by the QNMs related to 
V%. However, the above statements are possible only when there is a clear separation of these two interesting time 
scales. As we see in the next section, if we have a clear separation of the length scales M and I, this implies that 
Tm ~ 0(M). Thus a clear separation of the two interesting time scales — and - is then related to a clear separation 
of a c ~ l/l and ctb ~ 1/M. We note that when such a clear separation of the scales does not exist (as for e.g in the 
nearly extreme SdS black hole studied in 0,0) then we do not expect two different sets of QNMs. 

Let us return to (|33|) . Since we have not yet specialised to the case where the length scales are separated, we do 
expect the general solutions of this equation for uj to be complicated. However, if we specialise to the case where 
M -C Z, and work in an approximation where we can neglect O(Mfl) terms and higher, we must see solutions u that 
split into two sets as argued above. We therefore proceed to study the various terms in lj3H|l and drop terms that can 
be neglected in this approximation. For M <C Z, otb ~ 1/M and a c ~ l/l. Solutions uj for the case of the two length 
scales being well-separated are of the form uj ~ 0(1/M) orw~ 0(1/1). 

We now look for solutions of the form uj ~ 0(1/M). Then, in (|33|) . the first term has a factor (XV — X'Y). From 
the detailed expression for this factor, it is clear that this factor is proportional to exp(iwZln2). Now, uj = c<jq + iT 
where V > 0. Furthermore, T ~ 0(1/M). Therefore, this factor (XY' - X'Y) is proportional to exp(— Z/M). On the 
other hand, the second term in 1-SMI) has a factor (XX' — X'X) which is proportional to exp(Z/M). Since M <C I, the 
first term can be neglected compared to the second term (Here, we have also used the asymptotic expansion of the 
hypergeometric function for large parameter). Therefore, in this approximation we must solve the equation 

XX' - X'X = . (34) 

Here again, we look at the detailed expressions for X' in (|24|l and X' . Using these detailed expressions and multiplying 
(|34H by M, we see that the first term in (I34|l XX' is 0(M/l) compared to the second term X'X. Neglecting the first 
term, we see that we must now look for zeroes of the second term which lead to frequencies uj ~ 0(1/M). These are 
given by the zeroes of X' , which occur (as can be seen on inspecting l|24[l when 1 + a/2 = — n and n > is an integer. 
Thus, the frequencies uj are : 



w = ^^T-\+^+l)} ■ (35) 

We note, that above, we do not consider zeroes of X as they lead to frequencies not of 0(1/ M) as assumed while 
arriving at l|34|) . We have nevertheless shown that in the M <C I approximation, there indeed exist solutions u> ~ 
0(1/ M) to 

We now consider the other possibility, i.e w ~ 0(1/1). Now we can no longer neglect the first term of Il3.'i[l compared 
to the second term, as we did before. However, we multiply the whole equation by M, and observe that part of the 
second term can still be neglected as in the previous case of uj ~ 0(1/M). More precisely, the second term has a factor 
(XX' — X'X) where the first term is 0(M/l) compared to the second, and can be dropped. We now concentrate on 
the approximate equation 

(xf-x'Y) + x'x— 5r^r — ttt = - ( 36 ) 

r(ai)r(oi) T(a 1 - a + i)r(6 1 - a + l) 

There are more 0(M/l) terms left in the above equation. However, we observe that the first term has a factor l/r(ai). 
Thus zeroes of the above equation where u ~ 0(1/1) could occur at poles of V(a\) which are also zeroes of X. We 
check for such solutions. The poles of T(ai) occur when T(a\) = —n' and ml > is an integer. The corresponding 
frequencies are 




< 37 > 

The zeroes of X occur for a\ + 1 = —In where n > is an integer. The corresponding frequencies are 



'§-7 +<(2n+|)]. (38) 
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From l|37[) and l|38|l . we see that the poles of T(ai) which are also zeroes of X occur when, n' is an odd integer. These 
frequencies u> ~ 0(1/1) are then the solutions of (|36[l and in the M <C I approximation, the solutions of (|33|) . 



MAXIMUM OF THE SDS POTENTIAL 



In the previous section, we argued for the presence of two sets of QNMs, one proportional to ctb and the other to a c , 
in the approximation M C i. We also computed these frequencies which are given by (|35|l and (|38|l . But it remains 
to evaluate V m , the maximum of the SdS potential Vl{t). The maximum is a solution of the equation dVL/dr = 0, 
which is a sextic equation in r. For high values of the angular mode number L, the potential is approximately 

2M r 2 L(L + 1) , . 

^(r) ~ (1-—--)^-^. (39) 

Then the maximum occurs at r m = 3M. Thus, for large L, we can take V m = Vl(3M). 

x ,1 9Af 2 2 2 L(L + 1) S 

Let us now look at small values of L. We first consider the case L = 0. For L = 0, numerically, with M = 1 and 
A = 10 -4 , we can estimate the maximum, which is r m = 2.67. For M = 2 and the same value of A, r rn — 5.33. This 
suggests that r m — 2.67 M. Analytically, the equation to be solved is dVo/dr = 0, and the equation is 

r 6 Mr 3 3Mr 8A/ 2 , s 

2 ¥-^-nr + ir = - (41) 

Now, let us assume that the solution r m ~ O(M) as suggested by our numerical results. Then, in the equation above, 
the first and second terms are 0((M/l) 2 ) or higher powers of (M/l) compared to the last two terms. Therefore, 
dropping these terms we see that r m = |M. This agrees with our numerical results. For the QNMs with L — 0, we 



can therefore set V m = Vo(%M) 



For L = 1, we can again solve the equation dVi/dr = 0, assuming as before that the solution r m ~ 0(M) and in 
the approximation M <C I- We find that r m = ^(3 + y/7o)M, i.e r m ~ 2.886Af. We expect r m = 2>M to be a good 
approximation to the maximum for higher values of L. 

In all the above cases, we can explicitly check that V m /al and V m /a 2 are greater than 1/4. Thus, it is the real 
part of the QNM frequencies that depends on V m . Substituting for the value of V m from (|40l) in the expressions for 
the QNMs l|35|) and l|38|l - we get the explicit expression for the QNM frequencies for L > 0. For L = 0, using l|42|) 
for V m , we obtain more precise values for the QNM frequencies. 



VI. DISCUSSION 



In the previous sections, we have done an analysis of a perturbation in the potential V and given a sketch of its 
behaviour in time. We have found the presence of two sets of QNMs which are relevant at two different time scales 
(as discussed in Section IV): 



ui = 06[J-^.-i+i(2n+|)] and (43) 



a; 



h 



^-7 + «(2»+§)]. (44) 
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V m is determined as discussed in the previous section. We have also given a qualitative picture for the times at 
which each of these sets of QNMs would be relevant for the temporal behaviour of the scalar field. 

We now address the question of how well this reflects the time behaviour of a perturbation in the SdS potential. 
Numerical studies have been done for scalar fields in the SdS spacetime |7|, |8( ■ We indeed reproduce the most striking 
feature of the studies, the presence of two sets of modes which are proportional to the surface gravities of the black 
hole and cosmological horizons relevant at intermediate and late times respectively. 

Now let us make more quantitative comparison with the numerical results. We begin with the intermediate times 
at which the first set of modes (|35|) are relevant. In the numerical work a scalar field wave function with a 
Gaussian profile is chosen. Its subsequent time evolution in the SdS spacetime is found by numerically integrating the 
Klein-Gordon equation. A choice of r& = 1 and r c = 2000 (such that we are in the regime M <C I) is taken. A graph 
of the field behaviour in an SdS spacetime versus similar behaviour in a Schwarzschild spacetime is plotted for L = 
and L = 1. It is seen that for early times, the two graphs follow each other closely and display QNM oscillations - 
they start to differ from each other only for times much greater than 

Let us compare our first set of modes (|35|) with the Schwarzschild QNMs which seem to describe the behaviour in 
SdS well at these times. To do so, we can use an analytical approximation for the Schwarzschild QNMs by Ferrari 
and Mashhoon ;9j which approximates the QNMs well when their imaginary parts are not too large. Then in this 

approximation, the real and imaginary {ujf m ) parts of the Schwarzschild QNMs are wjj = a(J — |) and 

uij m — a(n+ i). Here, a and Uq are given by the height and curvature of the Schwarzschild potential at its maximum 
and n is an integer, a -1 = 3\/3-M. Comparing with the modes we see that for M <C I, Uq ~ V m . Also a^ 1 ~ AM 
and therefore ctb ~ (1.3)a. Thus we see that the real parts of the Schwarzschild QNMs and (|35|l have approximately 
the same M dependence. The angular mode L dependence is also the same in both real parts. However, coming to 
the imaginary parts, we see that our set of QNMs reproduces only those Schwarzschild QNMs starting from n = 2 
and subsequently reproduces only alternate n modes, the next one being n = 4. This probably has to do with our 
matching conditions at r m . 

We proceed to compare the time behaviour of the field in this model with numerical results for later times. In 
the next sub-section, we address intermediate time power law decay. The following sub-section compares the QNM 
behaviour at late times, i.e given by l|38|) in our model with numerical results - particularly with reference to angular 
mode dependence. 



A. Incompleteness of QNMs 



An important point in the numerical results 0, l£| is the observation of a power law behaviour in time for the field 
at certain intermediate times, rather than QNM behaviour. It is seen that the time decay of the field first displays 
Schwarzschild QNM behaviour, then Schwarzschild power law decay. At later times, the power law decay becomes 
faster (to be followed eventually by exponential decay). 

We examine whether there is a possibility of such interesting behaviour in our model. A detailed analysis has 
been done in (2l| on the conditions (and types of potentials) for which both QNM and power law behaviour in time 
may be present at different times. Also, the more general question of the potentials for which the QNMs form a 
complete set has been addressed. The results of the analysis are best expressed in terms of the Fourier transform 
of the Green's function for the system. This can be compactly written in terms of two auxiliary functions f(uj,x) 
and g(uj,x) which are solutions to the homogenous time-independent Klein Gordon equation, where / satisfies the 
QNM boundary condition at one boundary (in our case, the black hole horizon) and g satisfies the QNM boundary 
condition at the other boundary (in our case, the cosmological horizon). Then the Green's function is 

G(x,y;u) = f(u,x)g(uj,y)/W(uj) for < x < y , 

= f(u),y)g(u,x)/W(u) for < y < x . (45) 

Here, the Wronskian W(u>) = g(u),x)f'(u>,x) — f(u>,x)g'(u>,x) is independent of x. The QNMs are given by the 
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zeroes of the Wronskian, i.e when the functions / and g are linearly dependent. Following |2lJ, the QNMs do not 
describe the intermediate or late time behaviour completely when the functions g and / are not analytic w.r.t lu. 
Rather, non-analyticity of these functions could lead to power law behaviour in time. / and g are analytic in u> if 
the potential for the system V has 'no tail' at each of the boundaries x = — oo and x = +00 in the following sense 
(condition for boundary at x = 00) : 

/ dx x\V(x)\ < 00 , / dx xe ax \V(x)\ < 00 for any a > . (46) 
Jo Jo 

There is a corresponding condition for the boundary at x = —00. If this condition is violated for some a > ao > 0, 
then / and g may not be analytic in u> for Imoj > ao- Let us apply these results to our problem. For our problem, 
the function f(uj,x) obeying the QNM boundary condition at the black hole horizon is xii an d for x — > 00, it is given 
by the expression for Xqnm in 123- The function g(w,x) obeying the QNM boundary condition at the cosmological 
horizon is 

9M= (Q 2Pl -QiP 2 ) ^ 

It can be easily checked that the Wronskian W(u) is 

W(u>) = -2ioj(AQi + BQ 2 ) . (48) 

QNMs are frequencies for which the Wronskian is zero, and as expected, are given by solutions to (|32[) which we 
discussed extensively in the previous sections. 

We are interested in the question of whether / and g are analytic in u). The potential V for our problem goes to 
zero as e~ 2a " x as x — > 00 and as e 2abX as x — > —00. Therefore, conditions 146|) are violated for a > a& (coming from 
the condition at the black hole horizon) and for a > a c (coming from the condition at the cosmological horizon). We 
expect non-analyticity for two sets of u> : Imw > otb and Imio > a c . 

Let us now examine / and g to see if this is indeed the case, f = Xi- Using the properties of hypergeometric 
functions for the specific values of parameters a, b and c for the problem, we can write / in terms of the associated 
Legendre polynomial as 

xi = K(i - 0]^r( c )(£ - e) (1 - a - b)/A P^-i)/2^ ' 20 • (49) 

Therefore / is not analytic for c = —n where n > is an integer. This occurs for 

u) = iatb(n + 1). From l|T7|) . the poles of g are given by the poles of Q1/Q2, and therefore occur when ci = — n. The 
corresponding frequencies u> — ia c (n+ 1). Thus non-analyticity of / and g are as expected from the violations of l|46|l . 
This non-analyticity implies from the results in [2l| that QNMs do not completely describe the time behaviour of the 
scalar field at intermediate times for the potential V. Since we have chosen our potential to match the SdS potential 
asymptotically, we expect the same result for the SdS potential. The numerical results for a field in the SdS potential 
suggest power law behaviour in time for certain intermediate times. It would be interesting to see if such a power law 
behaviour could be derived in our model and the parameters the exponent would depend on. In particular, it would 
be useful if in our model, we could completely obtain the time dependence of the field and the times at which the 
power law behaviour sets in. This would be a valuable pointer to the time decay of the field in the SdS spacetime. 



B. Angular mode dependence of late time behaviour 



We now compare the late time behaviour in our model with numerical results. The late time behaviour in the SdS 
spacetime was first studied in p| and later in [8J. The numerical results consider (as for the early times) a scalar 
field wave function with certain initial conditions whose time evolution is determined by numerically integrating the 
Klein-Gordon equation. However, as it is required to access late times, it is no longer possible to work with the choice 
of black hole parameters used for early times, and the authors in [j| use r& = 1 and r c = 100 - but now the two 
horizons are not so well-separated. They then study the leading late time behaviour for the angular modes L = 0, 1, 2. 
For L = 0, they find that the field settles to a constant value. For L = 1,2, the field decays exponentially with 
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the argument of the exponent consistent with the formula —a c L. However, it is not possible for them to verify this 
proposed formula for the higher angular modes (L > 2) as the numerical integration becomes noisy. Thus the late 
time behaviour (being exponential decay) is a typical QNM behaviour where the QNM is pure imaginary. However, 
the numerical results display only the leading behaviour, given by the lowest QNM which seems to be cj = ia c L. It 
must be noted that no numerical results are available on the nature of the higher order modes for a given L (although 
the problem is addressed by the authors of in an analytical approximation where the central black hole in the SdS 
spacetime is disregarded). 

We consider the QNMs describing late time behaviour in our model given by (|38|l . These QNMs are proportional 
to a c and have a non-zero real part which depends on L. The dependence of the QNMs on the angular mode L is 
through the maximum of the potential V m and therefore affects only their real part. This is very similar to results 
obtained for the nearly extreme SdS black hole in 0, |(| . However it is at variance to the numerical results in Q 
where at least for the lowest QNM, the real part is zero, and the imaginary part depends on L (although both the 
model QNMs and the numerical results show a dependence on a c ). Unfortunately, we also have no numerical results 
available on higher order modes (for a given L) to compare with the higher order model QNMs. 

Our approximate potential V is the simplest model which is solvable in the limit of separation of scales, and reveals 
the presence of two sets of QNMs coming into play at intermediate and late times as seen in numerical studies. 
However, it does need to be refined further to also reproduce the L dependence of the second set of modes. However, 
the refinement has to be such that the model remains tractable at least in the physical regime of interest. This is a 
severe restriction. This is because much of our success above in reproducing the striking qualitative features has to 
do with correctly reproducing the asymptotic behaviour of the SdS potential at both horizons. The SdS potential 
falls off exponentially at both horizons, but with different arguments. This is difficult to reproduce in an exactly 
solvable potential, and therefore, one must piece together two potentials as we have done at a point, which could be 
the maximum of the potential. However, in this case, matching conditions may make the actual equation for QNMs 
complicated. The L dependence seen in numerical studies for the second set of modes and the failure of this model to 
reproduce them is related to the failure of a Poschl- Teller potential to describe the L dependence of the vacuum de 
Sitter potential. One obvious modification we could do to our model to reproduce the L dependence correctly is to use 
the same Poschl- Teller potential for r\, < r < r m and the vacuum de Sitter potential for r m < r < r c . However, in that 
case, the matching conditions at r m become very complicated. This makes the equation for QNMs complicated, and it 
is no longer possible to see the expected qualitative features emerge neatly in the M <C I approximation. Refinement 
of this model should therefore be concentrated on finding a tractable potential that is a good approximation to the 
vacuum de Sitter potential and could be pieced to the Poschl- Teller potential at the maximum. It should be however 
remembered here that piecing at the maximum, and resulting matching conditions may result in absence of certain 
QNMs from the spectrum. 



C. Further investigation of the model 



We have outlined above the steps to be taken for making the QNMs of this model agree well with numerical 
studies for late time behaviour of the field. However, it must be noted that no analytical approximation proposed 
so far for time behaviour of fields even in the Schwarzschild black hole reproduces its entire time behaviour from 
intermediate to late times. The analytical model in reproduces the QNMs, but not the characteristic late time 
tail for a Schwarzschild perturbation. It is to be noted that our model shows exponential late time decay for the field 
as seen for the SdS black hole, although the exponent does not match numerical studies. It also shows the QNMs at 
intermediate times. Therefore it offers an interesting platform to address the following questions: 

1) Can one map the complete time behaviour of a perturbation in this model, starting from intermediate times? 
The incompleteness of the QNMs of this model suggests interesting temporal behaviour when the field switches from 
a decay governed by the first QNM set to the second. It would be most exciting if in this model, a power law decay 
were found at these intermediate times. It would then shed light on similar behaviour in the SdS black hole. 

2) What happens to the near-extremal limit in this model? The discrepancies of our model QNMs with numerical 
studies are mainly in the angular mode dependence. We therefore expect that the results of taking the extremal limit 
in our model (two horizons approaching each other) should be a good approximation to the perturbation behaviour 
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when the SdS black hole approaches the Nariai black hole. Investigations of the near-extremal limit for asymptotically 
flat black holes (Reissner-Nordstrom) show pecularities |23 • Our model could be investigated for such pecularities in 
the near-extremal limit of the SdS black hole. 

It is possible to address the above two questions as this model appears tractable. For reasons elaborated above, 
such an investigation would indeed be valuable in understanding the response of the SdS black hole to perturbations. 
This investigation, and further refinements of this model to match numerical studies are the subjects of future work. 
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